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Введение. Во всех видах цифровой связи применяются методы помехоустойчивого кодирования. Во многих 
стандартах цифровой связи, например вай-фай (англ. \1-Е1) и 5@, используются коды с низкой плотностью 
проверок на четность. Эти коды популярны потому, что для них возможно построение кодеров и декодеров с 
невысокой вычислительной сложностью. Цель настоящей работы — повышение корректирующей способности 
известного битфлиппинг-декодера (англ. Б\-Ирриш®, ВЕ) ГОРС-кодов. Для этого строится модификация 
декодера, позволяющая динамически управлять одним из его основных параметров, выбор которого 
существенно влияет на качество декодирования. 

Материалы и методы. Рассмотрен известный декодер БИ-Иррше двоичных ГОРС-кодов. Некоторые его 
параметры не имеют жесткой связи с параметрами кода. С помощью имитационного моделирования 
исследована зависимость качества декодирования от выбора выходных параметров декодера Б\-Н1рришэ. 
Показано, что на результаты декодирования в этом случае существенно влияет входной параметр декодера — 
порог Т. Разработана модификация ВЕ-декодера двоичных ГОРС-кодов, в которой предлагается задавать порог 
динамически во время выполнения алгоритма в зависимости от степени повреждения кодового слова 
ошибками. Проведен сравнительный анализ корректирующей способности декодеров методом имитационного 
моделирования. 

Результаты исследования. Сформулирована и доказана лемма о максимальном значении порога Т декодера. 
Найдены верхние оценки для количества операций оригинального и модифицированного декодеров. Построена 
имитационная модель, реализующая цифровой помехоустойчивый канал связи. В модели исходные данные 
кодируются заданным ГОРС-кодом, зашумляются аддитивными равномерно распределенными ошибками, а 
затем поочередно декодируются алгоритмом Би-Е@ррше с различными параметрами порога Т и 
модифицированным декодером. По входным и выходным данным оценивается корректирующая способность 
использованных декодеров. Эксперименты показали, что в диапазоне реального уровня ошибок 
корректирующая способность модифицированного декодера выше, чем у оригинального, вне зависимости от 
выбора его параметров. 

Обсуждение и заключения. Доказанная в работе лемма устанавливает верхнюю границу значения порога в 
оригинальном декодере, что облегчает его настройку. По сравнению с оригинальным декодером разработанная 
модификация способна лучше исправлять ошибки. При этом сложность модификации увеличена незначительно 
по сравнению с оригинальным алгоритмом. Отмечено, что качество декодирования модифицированным 
декодером растет при увеличении длины кода и уменьшении количества циклов в графе Таннера, 
соответствующего проверочной матрице кода. 


Ключевые слова: ГОРС-коды, корректирующая способность декодера, динамический порог, двоичный 
симметричный канал, экспериментальное исследование. 
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Введение. В 1963 году в работе [1] Р. Галлагер впервые описал класс линейных блочных кодов, 
проверочная матрица которых содержит малое количество ненулевых элементов. Такие коды принято называть 
кодами с низкой плотностью проверок на четность или ГОРС-кодами (от англ. |10о\/-4епзИу рагцу сБесКк со4ез). 
Для них возможно построение кодеров и декодеров с невысокой вычислительной сложностью. Таким образом, 
при использовании ГОРС-кодов скорость передачи данных существенно не ограничивается. Многие 
современные работы посвящены ГОРС-кодам и их декодерам [2—5]. ГОРС-коды активно используются в 
разных стандартах цифровой связи, например вай-фай (англ. \/-Е!), 5С и оптической связи [6, 7]. Однако, 
несмотря на популярность этих кодов, некоторые связанные с ними задачи требуют исследования и решения. 
Одна из них — построение новых и улучшение существующих декодеров. 

Цель данной работы — повышение корректирующей способности известного декодера БИ-Иррше 
ГОРС-кодов (далее ВЕ-декодер). Для этого строится модификация декодера, позволяющая динамически 
управлять одним из его основных параметров, выбор которого существенно влияет на качество декодирования. 

Материалы и методы. Основными параметрами двоичных ГОРС-кодов являются длина М, 
размерность К и минимальное расстояние кода 4. Информационные слова [М№, К, 4]-кода С — это векторы 
т = (т, ть, ....тк) Е ЕХ, где Ё› — поле Галуа мощности 2, а кодовые слова — векторы с = (с1, со, ..., См) 
ЕЕ” [8]. Удобно задавать ГОРС-коды проверочной (М -К)жМ матрицей Н. Большее количество ее 
элементов — нулевые [1], поэтому удобнее хранить ее не целиком, а запоминая только позиции ненулевых 
элементов по строкам. 

Различают регулярные [9] и нерегулярные [10] ГОРС-коды. В регулярных кодах все строки и столбцы 
проверочных матриц содержат фиксированное количество единичных элементов (К и ] соответственно), иначе 
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код называют нерегулярным. Для удобства проверочные матрицы регулярных ГОРС-кодов будем называть 
регулярными матрицами, а нерегулярных ГОРС-кодов — нерегулярными. 

Регулярные ГОРС-коды обладают рядом преимуществ: легко оцениваемые параметры кода, удобство 
хранения матриц, невысокая вычислительная сложность алгоритмов кодирования и декодирования и пр. Кроме 
этого, декодеры регулярных кодов исправляют ошибки равномерно, в отличие от нерегулярных, которые 
исправляют ошибки в некоторых частях кодового слова хуже, чем в других. Однако задача генерации 
регулярных матриц с заданными свойствами является сложной, часто для ее решения используются методы 
перебора. 

Для обсуждения свойств матрицы Н удобно использовать соответствующий ей граф Таннера @ = 
(У, Е), где Е — множество ребер, а И = 5 УВ — множество вершин, 5 — множество строк матрицы Н, а В — 
множество ее столбцов [11]. Каждый ненулевой элемент Н, стоящий в 1-й строке и /-м столбце, задает ребро, 
соединяющее -ю вершину множества 5 и ]-ю вершину множества А. На рис. | представлен пример регулярной 


проверочной матрицы 3х6 с параметрами К = 4 и] = 2 и соответствующий ей граф Таннера. 
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Рис. 1. Цикл в графе Таннера и в проверочной матрице 


Верхний ряд вершин графа соответствует столбцам матрицы Н, а нижний ряд связан со строками Н. 
Важная характеристика проверочной матрицы Н ГОРС-кода — наличие и тип циклов в соответствующем ей 
графе Таннера. Цикл — это такая последовательность смежных вершин графа, в которой первая и последняя 
вершины совпадают. Длина этой последовательности называется длиной цикла. Минимальная длина цикла в 
графе называется обхватом. Если граф не содержит циклов, его обхват полагают бесконечным. Пример цикла 
длины 4 выделен жирными линиями на графе (рис. 1). 

Возможности по исправлению ошибок зависят не только от основных параметров ГОРС-кодов, но и от 
структуры проверочной матрицы Н. С одной стороны, наличие циклов небольших длин (таких как 4 и 6) 
заметно ухудшает корректирующую способность декодера [12]. С другой стороны, код, которому соответствует 
граф Таннера без циклов, не исправляет ошибок, т. к. его минимальное кодовое расстояние равно 2. Таким 
образом, задача построения проверочных матриц регулярных ГОРС-кодов является многопараметрической. 
При ее решении необходимо следить за основными параметрами кода, а также за циклами в графе Таннера, 
соответствующем проверочной матрице. 

Рассмотрим в удобном виде известный ВЕ-декодер ГОРС-кода С [13]. 


Вход: ГОРС-код С с параметрами [М№, К, а], заданный проверочной матрицей 


Й11 Й2 нии ам 
й й Бе й 
И 2 22 " м | (1) 
П(м-ю1 Им-к? + Им-Юм 
Вектор с’ =сС+е, СЕ С(С Е”), е(Е ЕЁ) — вектор ошибок; р — количество итераций алгоритма; Т — 


пороговое значение. 
Выход: кодовый вектор Сб Е С(< ЕЛ). 
Шаг 1. Положим счетчик 7 равным нулю. 
Шаг 2. Вычислим синдром 5 = С'НТ. Если $ = 0 илит = р, то переходим на шаг 5. 
Шаг 3. Выделим из вектора 5 = (51, 52,....5м_к) единичные координаты, т. е. $; = 1, = 1,(М-К). 


Составим множество [, = {1|5$; = 1}. Вычислим Й' = (11, А», ..., Ам), где 
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п; = ие, Ви. (2) 

Величины Ви, [= 1,..., М следует полагать неотрицательными целыми числами. Таким образом, А’ Е МХ, где 
№ = № {0}. 

Шаг 4. В векторе й' = (11,1», ..., Вл) находим все элементы й: > Т. Среди них выбираем случайный й! 
и инвертируем бит с! вектора с. Добавляем к счетчику 7` единицу и переходим на шаг 2. 

Шаг 5. С := С’. 

Исследования, проведенные в рамках данной работы, позволяют сделать ряд замечаний по ВЕ- 
декодеру. 

Замечание 1. Входной параметр р задает максимальное количество итераций алгоритма со 2-го по 4-Й 
шаги, но декодер может восстановить кодовое слово за меньшее число итераций. 

Замечание 2. При выборе параметра Т нужно руководствоваться следующими соображениями. Если 
известен параметр 4 используемого [М№, К, а]-кода С, то по нему можно вычислить Ё — число гарантированно 


исправляемых ошибок, и тогда количество итераций декодера ограничивается этим значением: 


а-1 
р=:= [> (3) 
Здесь |х| — округление числа х до меньшего целого. Если параметр неизвестен, то его можно оценить с 
помощью границы Синглтона [5] 
а<М-К+1 
и, используя (3), получить 
Ш к 
| 


Замечание 3. Структура декодера такова, что восстановление корректного кодового слова не 
гарантируется, даже если в зашумленном слове С’ = с + е возникло не более & ошибок (3). 


Замечание 4. В литературе для регулярных проверочных матриц в ВЕ-декодере рекомендуется 
] 
выбирать порог Т зависящим от веса ] столбца матрицы Н, а именно Т = >. Для нерегулярных матриц такие 


рекомендации в литературе не приведены. Корректирующую способность ВЕ-декодера может ухудшить 
неудачный выбор порога Т. При его большом значении на шаге 4 декодера в векторе А’ может не найтись 
координаты, превосходящей порог Т, следовательно, не будут исправлены ошибочные биты. При выборе 
малого значения Т на шаге 4 ВЕ-декодера в векторе А’ может появиться несколько координат, значение 
которых превышает порог. Среди них могут быть и координаты, не содержащие ошибку. Таким образом, выбор 
параметра Т может в значительной мере повлиять на качество декодирования. 

Результаты исследования. Сформулируем и докажем лемму о максимально возможном значении 
порога Т. Затем модифицируем ВЕ-декодер таким образом, чтобы порог устанавливался в процессе 
декодирования динамически, и проведем сравнительный анализ оригинального и модифицированного 
алгоритмов декодирования. 

Лемма. Пусть двоичный [М№, К, а]-код С задан с помощью проверочной матрицы Н, имеющей 
фиксированное количество / единичных элементов в каждом столбце. Тогда максимальное значение порога Т 
для ВЕ-декодера такого ГОРС-кода С не может быть больше 

Т=]- 1. (4) 

Доказательство. Пусть из канала передачи получен вектор С’ = С -е, где СЕС — верное кодовое 
слово, е Е ЕМ — вектор ошибок с весом Хэмминга и’(2). Если и’(г) = 0, то на шаге 2 вектор-синдром $ = 0. 
Следовательно, алгоритм перейдет на шаг 5 и вернет С’ в качестве ответа. В этом случае значение порога не 
используется. Если \’(е) > 0, то из регулярности Н вытекает справедливость неравенства Й/ <}, где й, — 
элементы вектора ЙА’. Инвертирование бит с, вектора С’ происходит в алгоритме, только если |! > Т. 
Следовательно, 

Т<# <]. 
Таким образом, формула (4) верна. 
Внесем в ВЕ-декодер изменения, которые позволят определять величину порога динамически, в 


зависимости от степени повреждения кодового вектора в канале передачи. 
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Вход: [М№, К, а]|-код С, заданный приведенной выше проверочной матрицей (1). Вектор с’ = С-+е, где 
СЕС(с Е), е (Е Е”) — вектор ошибок; р — количество итераций алгоритма; Т — некоторое пороговое 
значение, выбранное заранее. 

Выход: кодовый вектор сб Е С(с ЕЛ). 

Шаг 1. Положим счетчик 7 равным нулю. 

Шаг 2. Вычислим синдром 5 = С'НТ. Если $ = (0,...,0) илит = р, то переходим на шаг 7. 

Шаг 3. Выделим из вектора $ = (51,52, ....5м_к) единичные координаты, т.е. $; = 1, [= 1,(М-К). 
Составим множество [, = {15; = 1}. Вычислим В" = (№1, №»,...,Вм), где Ву такой же, как и в оригинальном 
декодере (2). При суммировании величины Йи [=1,...,М№М, следует полагать неотрицательными целыми 
числами. Таким образом, В’ Е №, где № = ММЦ {0}. 

Шаг 4. Инициализируем значение порога Т := тах(й1)1-1..м- 1. 

Шаг 5. Если Т > 0 

Выберем произвольный элемент Йа вектора р’ — такой, что №. >Т. 
Инвертируем бит са. 


Шаг 6. Добавим к счетчику Т единицу и перейдем на шаг 2. 

Шаг 7. с:=С.. 

Замечание 5. Модифицированный алгоритм в целом выполняет меньше итераций, чем ВЕ-декодер, т. к. 
на шаге 4 порог выбирается динамически. Следовательно, декодер не выполняет бесполезные итерации, на 
которых не изменяются биты вектора С’. Значение порога в модифицированном декодере зависит от числа 
ошибок, повредивших кодовое слово, и сразу устанавливается таким, что в зашумленное кодовое слово С’ 
гарантированно вносятся изменения. 

Оценим сверху количество операций сложения, сравнения и присваивания в обоих декодерах. В 
оригинальном ВЕ-декодере [М№, К, 4]-кода С производится р(КК + (М — КЮ)М + 1) операций сложения, р(3ЗМ — 
2К +2) операций сравнения и р((М — КЕ +2М+ 3) +1 операций присваивания. В ВЕ-декодере с 
динамическим порогом — р(КК + (М — К)М +3) операций сложения, р(5М№М — 2К + 3) операций сравнения и 
р((М —К)(К +1) +2М + 4) + 1 операций присваивания. Здесь р — параметр декодера, устанавливающий 
максимальное количество операций, К — вес строк проверочной матрицы кода. Заметим, что при реализации 
алгоритма фактически не используются операции умножения и деления, т. к. на втором шаге для вычисления 
синдрома $ удобно использовать операции сложения вместо умножения. Напомним, что матрица Н обладает 
разреженной структурой, и ее строки удобно хранить в виде списка номеров ненулевых элементов. 
Следовательно, вместо умножения вектора С’ на матрицу Н необходимо суммировать координаты вектора С’, 
номера которых совпадают с номерами ненулевых элементов в соответствующей строке матрицы Н. 

В сравнении с оригинальным алгоритмом модифицированный ВЕ-декодер выполняет больше 
операций, но не значительно: число операций сравнения увеличилось на р(2М + 1), присваивания — на р, 
сложения — на 2р. 

Для сравнительного исследования корректирующей способности оригинального и модифицированного 
алгоритмов декодирования создано программное средство, реализующее имитационную модель двоичного 
симметричного идеально синхронизированного помехоустойчивого канала связи согласно [14—16]. Для 
обеспечения помехоустойчивости в модели использованы ГОРС-коды и ВЕ-декодеры (оригинальный и с 
динамическим порогом). Ошибки в канале моделировались независимыми и равномерно распределенными. 

В экспериментах использованы специально найденные проверочные матрицы, задающие ГОРС-коды. 

Опишем основные параметры этих матриц, используя стандартные обозначения основных параметров кода, а 
также: ] и К — вес каждого столбца и вес каждой строки проверочной матрицы соответственно; 4, 6 —4иб 
циклов в графе Таннера, соответствующем проверочной матрице. 
Регулярная матрица Ну: М = 20, К = 5, =3, К = 4, 4 = 6, в. = 0, «в = 41. 
Регулярная матрица Но: М = 28, К = 7, ] =3, К = 4, а = 6, в. = 0, 5 = 42. 
Регулярная матрица Нз: М = 28, К = 7, ] =3, К =4, а = 6, в. = 0, вв = 29. 
Нерегулярная матрица Н»: М = 32, К = 5, ] = 3, а = 12, в. = 0, ®в = 0. 

С использованием этих матриц были построены кодеки ГОРС-кодов и проведены имитационные 
эксперименты. На рис. 2—5 представлены графики зависимости корректирующей способности построенных 
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кодеков ГОРС-кодов от вероятности ошибки в канале. Обоснование выбора значений порога Т =ТиТ=2в 


ВЕ-декодере см. в замечаниях 3, 4 и лемме. 
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Рис. 2. График корректирующей способности декодеровдля [ГОРС-кодов, заданных матрицей Н> 
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Рис. 3. График корректирующей способности декодеров для 1/)РС-кодов, заданных матрицей Нз 
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Рис. 4. График корректирующей способности декодеров для ГОРС-кодов, заданных матрицей Н. 
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Рис. 5. График корректирующей способности декодеров для ГОРС-кодов, заданных матрицей Н1 


В диапазоне реального уровня ошибок [8, 13, 14] на рис. 2—4 можно наблюдать, что ВЕ-декодер при 


Ь&р://уезбиК-доп$ва.га 


значении порога Т = 2 показывает лучшие результаты, чем при Т = 1, а модифицированный ВЕ-декодер 
обладает лучшей корректирующей способностью по сравнению с оригинальным. 

Декодеры демонстрируют схожую эффективность при малых значениях длины кода, однако при ее 
102 увеличении модифицированный декодер показывает лучшие результаты. Так, при вероятности ошибки в 
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непомехоустойчивом канале 0,05 различие в вероятности ошибки в помехоустойчивом канале при 
использовании ВЕ-декодера с порогом Т = 2 и Т=1 составляет от 0,005 до 0,03 в пользу использования 
большего значения порога. Если же задействовать ВЕ-декодер с порогом Т = 2 и модифицированный декодер, 
это различие колеблется в зависимости от [ГОРС-кода в интервале от 0,001 до 0,003. При вероятности ошибки в 
непомехоустойчивом канале 0,1 вероятность ошибки в помехоустойчивом канале при использовании ВЕ- 
декодера с порогом Т = 2 меньше, чем с порогом Т = 1, на величину от 0,001 до 0,02. При использовании ВЕ- 
декодера с порогом Т = 2 и модифицированного декодера это различие колеблется в зависимости от ГОРС- 
кода в интервале от 0,002 до 0,01. 

Оба декодера чувствительны к количеству циклов в графе Таннера, соответствующему проверочной 
матрице ГОРС-кода. Чем больше отношение количества циклов к общему количеству элементов в матрице, тем 
хуже исправляет ошибки любой ВЕ-декодер. При проведении экспериментов было интересно выяснить, можно 
ли так увеличить количество циклов в матрице, что модифицированный декодер покажет худшие результаты по 
сравнению с ВЕ-декодером. Опытным путем была найдена такая матрица — Н:, содержащая 41 цикл длины 6. 
Результаты исследования корректирующей способности декодеров для этой матрицы показаны на рис. 5. 
Заметим, однако, что матрица Н> содержит еще больше циклов длины 6, а именно 42. Принципиальное отличие 
матриц Н1 и Н> — в плотности единиц: 

— в Н1 — 60 единичных элементов на 300 элементов матрицы, 

— в Н> — 84 единицы на 588 элементов матрицы. 

Напомним, что особенностью ГОРС-кодов является разреженная структура проверочной матрицы, 
поэтому Н> более типична для ГОРС-кодов. 

Обсуждение и заключения. В работе рассмотрен декодер 6-Е ррше двоичных ГОРС-кодов. Даны 
рекомендации о выборе таких входных параметров декодера, как порог и количество итераций алгоритма. 
Сформулирована и доказана лемма о максимальном значении порога декодера. Разработана модификация ВЕ- 
декодера двоичных [ОРС-кодов, в которой предлагается задавать порог динамически во время выполнения 
алгоритма в зависимости от полученного синдрома. Для оригинального и модифицированного декодеров 
найдены верхние оценки количества операций. Эти оценки показывают, что модификация усложняет декодер 
незначительно. Проведенные имитационные эксперименты демонстрируют лучшую корректирующую 
способность модифицированного декодера по отношению к оригинальному. Эксперименты также показали 
зависимость качества декодирования от степени разреженности матрицы и количества циклов длины 6 в графе 
Таннера, соответствующего проверочной матрице ГОРС-кода. Таким образом, возникает задача построения 
проверочных матриц с малым количеством коротких циклов, что является предметом дальнейших 
исследований. 
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